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Abstract. The PHENIX experiment has shown that J/i/is are suppressed in central Au+Au collisions at a 
center of mass energy per nucleon-nucleon collision y^ijvjv = 200 GeV, and that the suppression is larger 
at forward than at mid-rapidity. Part of this difference may be explained by cold nuclear matter effects 
but the most central collisions suggest that regeneration mechanisms could be at play. In 2007, PHENIX 
collected almost four times more Au+Au collisions at this energy than used for previous published results. 
Moreover, the addition of a new reaction plane detector allows a much better analysis of the J/i/; behavior 
in the azimuthal plane. Since a large elliptic flow has been measured for open charm, measuring 3/^ 
azimuthal anisotropics may give a hint if J/i/j are recombined in the expanding matter. First PHENIX 
results of J ftp elliptic flow as a function of transverse momentum at forward rapidity are presented in this 
article. The analysis is detailed and results are compared to mid-rapidity PHENIX preliminary results as 
well as to predictions. 



PACS. 25.75.-q Relativistic heavy-ion collisions - 12.38.Mh Quark-giuon plasma 
flow - 13.20.Gd Decays of J/'i/', T, and other quarkonia 



25.75.Ld Collective 



1 Introduction 

A non-central collision creates a spatial anisotropy of the 
overlapping interaction region. In a thermodynamical pic- 
ture, the asymmetric distribution of the initial energy den- 
sity causes a larger pressure gradient in the shortest di- 
rection of the ellipsoidal medium which can lead to an 
anisotropic azimuthal emission of particles. Such an el- 
liptic flow is quantified by measuring the second Fourier 
coefficient, V2, of the produced particles in the (j) direction, 
with respect to the reaction plane angle, ip [Ull]: 
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A positive elliptic flow has been measured for light mesons [J 
which supports a picture with a rapid thermalization of 
the system. 

The cc pair that eventually forms a J/'0 is produced 
in the early part of the collision (through giuon fusion at 
RHIC energies). These are called direct J ftps. Part of the 
J/tfj yield also comes from feed down of excited states (■0' 
and Xc) and can contribute to up to ~ 40% jj. In heavy 
ion collisions, charmonia production rate can be reduced, 
for example, by gluon shadowing due to nuclear parton 
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distribution modification [5l[6], nuclear absorption [7| or 
color screening in a quark gluon plasma [5]. In addition, 
interactions with comovers [HI might also dissociate 3 ftps. 
On the other hand, given the large density of uncorrelated 
charm quarks in the hot medium at RHIC, these may 
recombine and enhance the 3/ip yield [Tnilll| . 

A way to experimentally disentangle direct from re- 
generated J ftps is to look for J/tp elliptic flow. Indeed 
PHENIX has measured a large elliptic flow for open charm 
from non-photonic electrons at mid-rapidity [12j as shown 
in Figure [T] for the same centrality class used for the J/ip 
V2 reported here. If regeneration is strong in this rapidity 
region, one would expect 3 ftps to inherit the positive el- 
liptic flow of the open charm from which they are formed. 
Assuming similar open charm elliptic flow at forward ra- 
ipidity, measuring 3/'ip V2 at various rapidities could in- 
dicate if the difference seen in their nuclear modification 
factor comes from regeneration. If no J/'0 elliptic flow is 
measured, there is little chance that the measured 3/ips 
come from regeneration. 



2 Detectors 

Elliptic flow is measured with respect to the collision reac- 
tion plane estimated using a subset of particles produced 
during the collision pj[T3] . Careful attention must be taken 
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Fig. 1. Non-photonic electrons V2 as a function of px for the 
centrality 20-60% at mid-rapidity {\y\ < 0.35) 



to choose the subset so that all correlations other than the 
elliptic flow itself are removed between this subset and the 
particle of interest. Until 2007 the collision reaction plane 
was measured in PHENIX by using Beam-Beam counters 
(BBC). In 2007 a new reaction plane detector (RxnP) was 
installed. It allows two times better precision on the reac- 
tion plane measurement, as shown in Figure [2l The res- 
olution by which the reaction plane angle is measured is 
used as a correction to the measured V2 to account for the 
finite precision of the reaction plane measurement. 
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Fig. 2. Reaction plane resolution correction {{cos{2A'4)rp))) 
as a function of centrality, measured with the RxnP in circles, 
the Muon Piston Calorimeters (MPC) in triangles, or the BBC 
in squares. Larger values correspond to higher precision in the 
reaction plane determination. 



At forward rapidity, i /ips are measured through their 
decay into di-muons with two spectrometers, one on each 
side of the interaction point. Each spectrometer is com- 
posed of muon identification (MuID) and tracking (MuTr) 



systems [M]. The MuID, made of larroci tubes and steel 
absorbers, allows the identification of the muons through 
their penetration depth. The MuTr is made of cathode 
strip chambers that measure the particles momentum through 
their bending in a magnetic field. The muon spectrometers 
have a rapidity coverage for single muons of 1 77 1 G [1.2—2.2] 
which overlaps with the RxnP acceptance (|?7| G [1 — 2.8]). 
To eliminate possible bias from radiated gluons accom- 
panying i/ips at similar rapidity, only the RxnP detector 
opposite to the arm where the muons go is used in this 
analysis. 



3 Analysis Method 

The integrated luminosity used for the analysis for this 
result is 537 /ib^^ (611 ^b~^) at negative (positive) rapid- 
ity, which is almost four times more than in the previous 
publications |15| . This analysis is performed on an online- 
filtered data sample which enabled results very soon after 
the data was taken. The online filtering was based on fast 
tracking in the MuID and has no effect on the V2 derived 
here. 

The elliptic flow is obtained by dividing J /■i/;s into only 
two (/> — ?/' bins: ]0,7r/4] and ]7r/4,7r/2]. V2 is then calculated 
by comparing the number of J /ips that belong to the bin 
which contains the reaction plane, A^'", with those that 
are in the bin which does not contain the reaction plane 
jyout^ following the formulae: 



(iVin _ 7v°"') 
{N"^ + iV°"') 



(2) 



with uncertainty. 



7r/2 



where cr™ is the error on A^'" and cr°"* is the error 
on A^°"'. Finally, v™'^^^ is divided by the reaction plane 
resolution in order to get the true J/t/j V2 as described 
in [13]. 

Since the background subtraction is critical for the 
V2 measurement, improvement on the subtraction method 
has been made for this analysis. Previously, a mixed-event 
subtraction technique was used to extract the J/ip signal 
in a high background environment |15| . This method ac- 
counts for the fact that the MuTr acceptance differs for 

like-sign (++or ) and unlike-sign (H — ) muon pairs and 

improves the statistical errors for bins where the signal 
over background is poor. The invariant mass distribution 
using the mixed-event subtraction for the 20-60% central- 
ity range is shown in Figure [31 The signal fit seems to be 
accurate, but the spectrum is distorted for masses below 
2.6 GeV/c^. The distortions originate from a bias intro- 
duced by the online filtering which creates a difference in 
the mixed-event (Mixed) and same-event (FG) distribu- 
tions. To account for this bias, we use both like-sign and 
mixed-event distributions as follow: 



Mixedj 
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3 



The like-sign distributions help correcting for acceptance 
effects caused by the online filtering and the mixed-event 
distributions allow to reproduce a statistics unlimited back- 
ground. 

The mass spectrum obtained with the combined sub- 
traction is shown in Figure ID The low mass distortions 
are highly reduced and the residual background below the 
J/V" peak should also be reduced. On the other hand, this 
method results in statistical errors that are larger by a 
factor of about V2 due to the use of the unlike-sign same- 
event distributions. 
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Fig. 3. 20-60% invariant mass distribution using the mixed- 
event background subtraction. The signal is obtained averaging 
results of three different fits |15) . 
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Fig. 4. 20-60% invariant mass distribution using the combined 
background subtraction. The signal is obtained averaging re- 
sults of three different fits [15j . 



4 Results 

The elliptic flow at forward rapidity {\y\ £ [1.2 — 2.2]) as a 
function of pT for the 20-60% centrality selection is shown 
in Figure [5] for negative (positive) rapidity in black circles 
(red squares). In Au+Au collisions, the forward-backward 



rapidity symmetry allows one to average the two mea- 
surements. The averaged result is shown by the magenta 
closed circles in Figure|6l The error bars on the two figures 
account for statistical and point-to-point uncorrelated er- 
rors. They come from the statistical uncertainty on the 
number of signal counts [15| . and the systematic uncer- 
tainty of the signal counting and background line shapes. 
The boxes represent point-to-point correlated errors that 
account for the error on the average RxnP detector resolu- 
tion and the error on the J/ip (j) angle measurement (least 
accurate for pT < 1 GeV/c). An additional global relative 
uncertainty of 3%, written on the figures, accounts for 
the error on the technique used to determine the reaction 
plane angle and resolution [2l[T3]. 

^ O.Si I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i_ 
„ ^ £^[0,5] south ^,= -0.066 ± 0.1 44 ±IZ • |V - -l^^l^ 

= northv,= -0.105±0.146±°°°= " = 

0.3+ + 

= Global relative syst. ± 3% I = 

0.2+ + 

0.1 + 1 + 

= =•= = 

-0.1 + ™p + 

= ^ m = 

-0-2^ ^ centrality [20,60%] ^ 

-0.3 + opposite RxnP ^ 

.(,4+ online filtered J/i|f^uV* + 

E PHEI^IX prqlimina|-Y , , , i , , , , i , , , , i , , , , i^^t^*^ ?P9,9^y ,= 
■"■^0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Pt 

Fig. 5. V2 as a function of pr for negative rapidity arm in 
black circles (red squares) at positive rapidity. 
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Fig. 6. V2 as a function of pt for combined forward rapidity in 
magenta closed circles (full data sample) and at mid-rapidity 
in open blue circles (42% of the data sample) , with theoretical 
predictions made for mid-rapidity. 

Also shown in Figure [6] are PHENIX preliminary mea- 
surements at mid-rapidity {\y\ < 0.35) (TB] in open circles. 
Forward and mid-rapidity measurements are independent 
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measurements since they use different detectors and trig- 
gers. In addition the methods to extract V2 are also dif- 
ferent and thus so are the uncertainties. At mid-rapidity, 
the point-to-point correlated systematic uncertainties are 
dominated by the unknown shape of the background V2- 
At forward rapidity, only corrections for the finite preci- 
sion of the reaction plane measurement and angle mat- 
ter, after the total subtraction of the background. Forward 
and mid-rapidity measurements are perfectly compatible 
with each other for all pT bin. The V2 obtained for each 
rapidity for pT e [0 - 5] GeV/c is: -0.094 ± 0.104±;] o;]f 
at |y| e [1.2 - 2.2] and -0.10 ± 0.10 ± 0.02 at |y| < 0.35. 

It is unclear whether i/tp elliptic flow at forward ra- 
pidity should differ from mid-rapidity, especially since the 
forward rapidity measurement only reaches \y\ = 2.2 and 
the underlying physics related to collective behavior might 
not differ much over such a rapidity range. Given this, an 
option is to combine the two rapidity measurements to in- 
crease statistical significance of the results. The statistical 
compatibility for V2 to be positive for G [0 — 5] GcV/c 
is only of the order of 10%. It is dominated by the [1 — 
2] GeV/c bin which is the bin where the J/t/'S pt distri- 
bution is peaked. The probability of having a positive V2 
in [1 — 2] GeV/c is only 6%. It is more interesting to know 
if V2 is positive and to what extent for pT > 1 GeV/c 
since coalescence models predict an increase of J/V' flow 
in this region as seen in Figure [6l The statistical proba- 
bility for the combined V2 measurements to be positive in 
[2 - 5] GeV/c is 42%, and for it to be above 0.1 is 15%, 
which does not allow to distinguish between models. Ad- 
ditional data especially for px > i GeV/c would give more 
statistical significance. 

Several predictions of V2 are shown in Figure [6] Higher 
V2 is predicted when more recombination is at plav[T7 ( fT8 l 
[inilinj, and almost no V2 for direct J/V'S POII^ . even when 
including absorption and comover interactions |22| . It is 
to be noted that these predictions are for minimum bias 
collisions, and computed for mid-rapidity. No predictions 
at forward rapidity are yet available. 

The J/ip V2 measured by PHENIX can be but into per- 
spective compared to recent lower energy measurements 
at the SPS. NA60 measured a positive V2 of 7 ± 3% 
in In+In non-central collisions with no selection in 
Pt- At the SPS energy and in In-|-In collisions, given the 
much smaller charm cross section, it is unlikely that a 
sufficient number of cc pairs is produced to allow substan- 
tial formation of J/V' by recombination. If confirmed, the 
measured positive flow at the SPS is likely caused by other 
phenomena. One might speculate that nuclear absorption 
is responsible, despite earlier expectations that this effect 
should be negligible. Also at RHIC energies there may be 
additional contributions affecting v2. In any case, models 
should provide predictions for different nuclei and ener- 
gies so that the measurements can distinguish any other 
effects from elliptic flow induced by recombination in a 
QGP. We note that NA50 measured a small 3/ip V2 in 
Pb+Pb collisions [24] with a maximum of 3.5 ± 1.5 ± 1.3% 
for Pt ^ 2 GeV/c. This result is not binned in centrality, 
and since V2 for central collisions should be close to zero. 



their measurement could correspond to slightly larger val- 
ues for non-central collisions. Interpretations of these re- 
sults would also help to shed light on higher energy mea- 
surements done at RHIC. 



5 Outlook 

PHENIX has measured J/i/; azimuthal anisotropy at for- 
ward rapidity for a centrality of 20-60% using all the statis- 
tics from the 2007 data taking, and an online-filtered sam- 
ple. Forward and mid-rapidity results are compatible within 
large uncertainties. For each measured pt bin, the value 
is compatible with both zero and maximal flow. The cur- 
rent precision on the measurements does not allow firm 
conclusions. 

The forward rapidity statistical significance will be im- 
proved for final results by using the minimum bias rather 
than the online-filtered sample, with a gain in statistics by 
about ~10% and by benefiting from the improved han- 
dling of the combinatorial background described earlier, 
leading to a reduction of the statistical uncertainty by ^/2. 
A much larger data sample is expected to be available at 
RHIC in 2010. 
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